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Gravitational parity violation is a possibility motivated by particle physics, string theory and 
loop quantum gravity. One effect of it is amplitude birefringence of gravitational waves, whereby 
left and right circularly-polarized waves propagate at the same speed but with different amplitude 
evolution. Here we propose a test of this effect through coincident observations of gravitational waves 
and short gamma-ray bursts from binary mergers involving neutron stars. Such gravitational waves 
are highly left or right circularly-polarized due to the geometry of the merger. Using localization 
information from the gamma-ray burst, ground-based gravitational wave detectors can measure the 
distance to the source with reasonable accuracy. An electromagnetic determination of the redshift 
from an afterglow or host galaxy yields an independent measure of this distance. Gravitational 
parity violation would manifest itself as a discrepancy between these two distance measurements. 
We exemplify such a test by considering one specific effective theory that leads to such gravitational 
parity-violation, Chern-Simons gravity. We show that the advanced LIGO- Virgo network and all- 
sky gamma-ray telescopes can be sensitive to the propagating sector of Chern-Simons gravitational 
parity violation to a level roughly two orders of magnitude better than current stationary constraints 
from the LAGEOS satellites. 
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I. INTRODUCTION 

Several reasons exist to believe that there is some de- 
gree of parity violation in gravity, at least at the quan- 
tum scale and possibly carrying through to macroscopic 
scales through the anomaly mechanism. From an experi- 
mental standpoint, it is curious that the weak interaction 
exhibits maximal parity violation, yet the other funda- 
mental forces apparently exhibit none. From a theoreti- 
cal standpoint, parity violation arises from the standard 
model of particle physics [l|-[3| , some sectors of string the- 
ory [4j, [5| , and extensions of loop quantum gravity |6l4l0| . 

Gravitational parity violation can be encoded in the 
action as an extension to general relativity (GR) that 
consists of the addition of an antisymmetric product 
of curvature tensors coupled to a scalar field, which to 
leading order in the curvature leads uniquely to Chern- 
Simons (CS) gravity [llj. Jackiw and Pi [12| demon- 
strated that the CS term can be covariantly embedded 
with its three dimensional counterpart as a consistent 
modification to GR. In the standard model, there ex- 
ists a radiatively generated one-loop chiral anomaly that 
gives rise to a CS term coupled to lepton number [13[. In 
heterotic and Type I supcrstring theories, the CS term 
is generated through the Green-Schwarz gauge anomaly- 
canceling mechanism [14j |. In loop quantum gravity, the 
scalarization of the Barbero-Immirzi parameter coupled 
to fermions has been shown to lead to an effective CS 
action p-[lO[ • Even devoid of a specific fundamental the- 
ory, one can also show that the CS term unavoidably 
arises in effective field theories of inflation as one of the 



few non- vanishing, second-order curvature corrections to 
the Einstein-Hilbert action [15[ . Another phenomenolog- 
ical parity-violating extension to GR has been proposed 
where the fundamental constants of nature (e.g., New- 
ton's constant) break parity invariance [16| . 

Most investigations of gravitational wave (GW) signa- 
tures of parity violation have focused on amplitude bire- 
fringence. That is, compared to GR, right circularly- 
polarized waves are enhanced or suppressed and left 
circularly-polarized waves are suppressed or enhanced 
as they propagate [0, HI, [l7H19J . Such birefringence 
occurs when GWs propagate on a flat or curved back- 
ground. In fact, the propagation of such waves over 
cosmological distances has been proposed as an expla- 
nation of the baryogenesis problem in the early Uni- 
verse [20| and could have observable effects on the cosmic 
microwave back grou nd fla \2VL [221 and a stochastic GW 
background [23 



24l | . Cosmological GW amplitude bire- 



fringence has been discussed in the context of GW obser- 
vations of cosmological supermassive binary black holes 
(BHs) with the planned Laser Interferometer Space An- 
tenna (LISA) [l8J|. 

Amplitude birefringence arises in propagating modes 
if the parity operator does not commute with the Hamil- 
tonian. In its most general form, the effect of the plane- 
wave propagator can be expressed as 
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Here uj is the GW angular frequency, t is some time co- 
ordinate, /i+, x ,fe are the Fourier components of the GW 



of wavcnumber k, and u and v are multiplicative fac- 
tors characterizing the translation invariant amplifica- 
tion/suppression of each Fourier mode during propaga- 
tion in time. The quantity u corrects for background cur- 
vature effects: u = 1 in a flat-background, but it acquires 
redshift corrections in a cosmological background (see 
e.g. j25(). Given a specific modified theory, Eq. ([T]) can be 
obtained by solving the modified field equations for the 
propagating modes of the metric perturbation, where v 
measures the degree of parity violation. The right- and 
left-circular polarizations h B L = (h + ± ih x )/y/2 then 
propagate as 
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In the presence of amplitude birefringence, the eigenvalue 
matrix of the propagator operator remains diagonal, but 
with modified eigenvalues. Clearly then, right circularly- 
polarized waves are amplified or suppressed, while left 
circularly-polarized waves are suppressed or enhanced 
relative to GR, depending on the sign of v. 

The amplitude birefringence discussed above and in 
the rest of this paper is a pure propagation effect, which 
accumulates over propagation distance, in contrast with 
wave generation effects which do not. As the latter cause 
parity violation "at the source," their effect is observable 
at any distance [Lg, [26| . Wave generation modifications 
also couple to matter as well as curvature [27] , for exam- 
ple leading to weakly composition dependent waveforms 
in BH-NS versus NS-NS mergers. Recently, HH stud- 
ied parity violation in wave generation for extreme-mass 
ratio inspirals, finding that indeed trajectories and the 
multipolar generation scheme itself encode such viola- 
tion at a fundamental level. The influence of such wave 
generation effects on GW obscrvablcs would require a 
full (non-extreme mass ratio) post-Newtonian and mul- 
tipolar analysis, which is beyond the scope of this paper. 
Instead, here we concentrate on wave propagation effects 
only and for simplicity adopt pure GR wave generation 
without modificd-gravity corrections, recognizing that (i) 
such wave generation effects are calculable on a thcory- 
by-theory basis and (ii) their distance-independent in- 
fluence on binary phase should be easily distinguishable 
from any secular amplification accumulated over cosmo- 
logical distances. 

In this paper, we show that an interesting test of grav- 
itational parity violation due to propagation effects can 
be made with existing instruments: the SWIFT [29| and 
GLAST/Fermi [30[ gamma-ray satellites, and ground- 
based LIGO [3l| and Virgo |32| GW detectors after on- 
going upgrades. Such a test relies on the detection of co- 
incident GW/gamma-ray burst (GRB) events, a possible 
progenitor of which are neutron star (NS) binary merg- 
ers, that lead to a redshift observation either from after- 
glow or host galaxy identification. The reasons why such 
events are excellent to constrain parity violation are the 
same as those outlined in work that proposes the use of 
GW/GRB events to constrain GRB parameters [33[ and 



cosmological parameters [25J, |3J, [3Jf . Due to the collima- 
tion of the jet producing the GRB, the binary's orbital 
angular momentum at merger must point along the line 
of sight to Earth, and thus the GW signal must be highly 
circularly-polarized. Together with the accurate sky lo- 
cation obtained from the GRB, this allows an accurate 
distance measurement from the GW signal. This dis- 
tance measurement can then be compared to the purely 
electromagnetic distance measurement, which for typical 
redshifts z ~ 0.1 depends only on the Hubble constant. 
Amplitude birefringence will manifest as a discrepancy 
between these two distance measurements. If no discrep- 
ancy is found, the error on the distance measurements 
can be used to place upper limits on the possible degree 
of gravitational parity violation. 

The test proposed here is unique in that it employs 
only the radiative sector of any generic, gravitational 
parity violating theory. Current constraints on gravita- 
tional parity violation concentrate on a specific alterna- 
tive theory (CS gravity) and explore only its stationary 
sector through Solar System experiments [36I438J and bi- 
nary pulsar tests |39| ■ As sucn ; these tests constrain only 
specific and local deviations from gravitational parity in- 
variance in the neighborhood of the gravitational source 
- that is, near Earth in the Solar System case and near 
J0737-3039 in the binary pulsar case. The test we discuss 
here constrains generic parity violation over distances of 
hundreds of Mpc, along the light-cone on which the am- 
plitude birefringent GWs propagate. 

This test also compares favorably to future GW tests 
with LISA [ljl LLSj]. The fractional error in a distance 
measurement is inversely proportional to the signal-to- 
noise ratio (SNR) p. Although the latter is larger in the 
LISA case, the sensitivity to the parity violating effect 
is greater here, because this effect also depends on the 
GW frequency, which is much larger during NS merg- 
ers. Propagation effects should accumulate with distance 
propagated, and on dimensional grounds this distance 
should be measured in gravitational wavelengths - that 
is, the effect should be proportional to Df, where D 
is the distance propagated and / is the GW frequency. 
While GW/GRB events will have smaller Ds and lower 
ps than LISA supermassive BH observations (both by 
1-2 orders of magnitude), / is 5-6 orders of magnitude 
higher for the LIGO- Virgo network. Thus, the sensitiv- 
ity to parity violation could be 2-3 orders of magnitude 
better than in the LISA case, although as we shall see 
later this is somewhat reduced by systematic errors. 

The straw-man chosen here to compare against ex- 
isting tests of gravitational parity violation is CS grav- 
ity. In this theory, the magnitude of gravitational par- 
ity violation is controlled by a length parameter (in ge- 
ometric units), proportional to certain time derivatives 
of the CS coupling field 9. Solar system observations 
of the LAGEOS satellites [38( have placed the bound 
6q < 2000 km at 2a when 9q = 0. (The subscripts refer to 
evaluation at present, i.e. zero redshift.) We have found 
that the above constraint can be generalized assuming 
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FIG. 1: Given a coincident right-handed GW/GRB detection 
at an SNR of 10, the gray-shaded region (red in the color ver- 
sion) shows the area of (8o,0o) parameter space that could 
be excluded at 2a by a measurement of one right-handed 
short GRB merger waveform consistent with its expected 
GR amplitude, given host galaxy distance and (measured) 
chirp mass. The region excluded by Solar System observa- 
tions (|0o - 0o/H o \ < 2000 km at 2a) is outside the dotted 
lines. 



9 evolves on cosmological timescales. Allowing for both 
(#o,#o) to be nonzero to leading order, we can general- 
ize the Solar System constraint via the replacement rule 
Gq — > 1 0o ~~ 9q/Hq\, which then leads to the much more 
powerful bound |0 O - O /-Ho| < 2000 km. 

Our calculations suggest that the LIGO- Virgo network 
will be dramatically more sensitive to parity violation, 
as we show in Fig. Q] for a right-handed GW observa- 
tion. The gray-shaded regions (red in the color version) 
show the areas of phase space that could be detected 
or excluded with a coincident LIGO GW/GRB obser- 
vation at 2cr. This is to be compared to the regions of 
phase space ruled out by the generalized Solar System 
constraints, denoted in the figure as the region outside 
dotted lines. Inside the dotted lines, the generalized So- 
lar System constraint cannot rule out the CS modifica- 
tion. A GW observation could place much tighter con- 
straints on the allowed region of parity violating phase 
space. Such constraints are interesting because gravita- 
tional parity violation need not be Planck suppressed in 
the fundamental theory from which CS gravity derives. 
This is because the effective CS theory is shift-invariant, 
and thus, the scalar field that sources the correction must 
be massless, if this symmetry is not to be broken. 

A certain asymmetry in the possible constraints 
we could place on (0o,0o) given a single, coincident 
GW/GRB detection is clearly discernible in Fig. [JJ This 
asymmetry is an artifact of assuming a purely right- 
polarized GW event. In CS gravity, right-polarized sig- 
nals possess a resonant behavior if 9q > for frequencies 
in the LIGO sensitivity band, which greatly enhances the 
amount of parity violation and the constraints one could 
place. This situation, however, reverses if instead we had 



assumed a left-polarized GW observation. In such a case, 
we would have obtained exclusion regions identical to 
those of Fig. [T] but reflected under (0o,0o) — > (— &o, — 0o)- 
Therefore, a two-detector observation of a right and a 
left-polarized GW would allow us to constrain (0o,0o) 
equally well for both positive and negative 6$, thus ex- 
tending the exclusion regions of Fig. [TJ 

All the results we just described assume a somewhat 
idealistic test, where we have ignored certain astrophys- 
ical uncertainties, such as inaccuracies in host galaxy 
identification. Perhaps one of the most important un- 
certainties relates to the inclination angle. GRB jets are 
not necessarily perfectly aligned with the line of sight, 
which implies that the inclination angle is not exactly 
known a priori and the GW signal is not perfectly circu- 
larly polarized. Since the inclination angle is degenerate 
with the luminosity distance in the GW response func- 
tion, this deteriorates GW measurements of the latter by 
roughly a factor of two. However, even after taking sev- 
eral such uncertainties into account, we find that 8q could 
be constrained at the 10-50 km level, which is still much 
better than current Solar System constraints and even 
LISA. Moreover, even when accounting for such uncer- 
tainties, coincident GW/GRB measurements could con- 
strain a much larger region of the parity violating phase 
space all along the light-cone in which GWs propagate. 

The remainder of this paper presents additional details 
that lead to the results described above. We begin with 
a description of CS gravity and amplitude birefringence 
in Sec. [TT] as a warm-up for the introduction of generic 
gravitational parity violation and its GW observable in 
Sec. IIIII We then proceed with a somewhat idealized 
test in Sec. IIV1 where we compute the modified SNR for 
a parity-violating signal. Assuming a perfect electromag- 
netic and GW determination of the luminosity distance, 
we then compute the constraint one could place on parity 
violation, assuming the latter can be treated as a linear 
correction to GR. In CS gravity, however, certain reso- 
nances arise in the GW signal at the non-linear level when 
parity violation is allowed to be strong. Such resonances 
greatly enhance the parity-violating corrections to the 
SNR, and thus, allow one to place the much stronger con- 
straints shown in Fig. [TJ as we show explicitly in Sec. [\H 
One might worry that resonances could signal the pres- 
ence of instabilities, but we show that this is not the case 
in the Appendix. The proposed test, however, is rather 
idealized because correlations with other parameters and 
astrophysical uncertainties generically degrade measure- 
ments of the luminosity distance. These potential sources 
of error are discussed in Sec. I VII and combined with the 
idealized analysis to construct realistic estimates of the 
possible constraints one could place on parity violation. 

In the rest of this paper we employ geometric units, 
where the speed of light and Newton's gravitational con- 
stant are set to unity c = 1 = G. We further follow the 
conventions of Misncr, Thorne and Wheeler j40|, where 
the metric signature is (—,+,+, +), spacctimc indices are 
denoted with Greek letters, spatial indices with Latin let- 



ters and abstract tensors are denoted in boldface type. 



II. CHERN-SIMONS MODIFIED GRAVITY 

A. CS Basics 

CS gravity postulates the addition of a parity- violating 
term to the vacuum Einstein-Hilbert action (see [11| for 
a recent review) such that [l2[ 



S = 
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where g is the determinant of the spacetime metric, R is 
the Ricci scalar, and 
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in terms of the Ricmann and Levi-Civita tensors R and 
e. Such a term is explicitly parity violating because after 
a spatial triad inversion, the Levi-Civita tensor changes 
sign, while the square of the Riemann does not. There- 
fore, the CS correction to the action, the second term 
in Eq. ([3]), also changes sign under parity inversion, pro- 
vided 9 is a scalar field, i.e., it is even under parity. 

The CS field equations add terms involving V# and 
V 2 6* to the Einstein field equations: 



G + C = 8ttT, 



(5) 



where G is the Einstein tensor, T is the stress energy 
tensor and C is the C-tensor 
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We consider here both canonical and dynamical CS 
gravity. In the case of the former, we let 9 be an arbitrary 
function of time, such that it possesses the same symme- 
tries as a Friedman-Roberston- Walker (FRW) spacetime 
(homogeneity and isotropy). In the case of the latter, the 
9 field must satisfy its equation of motion, i.e., Eq. ([?])■ 
whose solution contains a homogeneous and an inhomo- 
geneous piece. The inhomogeneous solution is sourced by 
R*R, which is an extremely small quantity of 0(h 2 ) [28[ 
in the wave- or radiation-zone (many gravitational wave- 
lengths away from the source), where h is the norm of 
the metric perturbation. Therefore, the inhomogeneous 
solution is also small at the location of the observer and, 
in fact, parity odd because it must possess the same sym- 
metries as R*R. Since such a parity odd 9 would lead to 
a parity preserving correction, as then the combination 
6>R*R is parity even, we neglect the inhomogeneous so- 
lution henceforth. On the other hand, the homogeneous 
solution forces 9 to be a freely propagating wave, but 
since we require 9 to be homogenous and isotropic, so 
that the background symmetries are satisfied, we shall 
assume it is a free-function of time only. 

At this junction, let us make a brief detour to discuss 
the the units used for the scalar field. As written abovc ; 
9 has units of squared length, and its time derivative 9 
defines a length scale (or an inverse energy or mass scale 
in natural units), which is proportional to the coupling 
strength of the CS term. To make contact with stan- 
dard notation [llfl . 9 is here equivalent to (167rG)9a 2 / f3 
in the dynamical case or (16irG)9a in the canonical or 
non-dynamical case, where (a, /3) are CS coupling con- 
stants (see eg. fill]). 



B. Gravitational Waves in CS gravity 



and parenthesis stand for index symmetrization. The 
modified field equations arc this simple because the CS 
correction can be thought of as a boundary term that 
acquires dynamics due to the scalar field coupling [121 ] . 

CS gravity can lead to two somewhat distinct theories, 
depending on whether the CS scalar field is allowed to 
evolve dynamically or not. In what we shall call "canon- 
ical" or "non-dynamical" CS gravity [12| 9 is prescribed 
externally, and a common choice ("static CS gravity") is 
to let V6* be a constant vector pointing in the direction 
of time defined by cosmic expansion. In what we shall 
call "dynamical" CS gravity [3, |38|, |41j 9 is promoted to a 
dynamical field by adding a kinetic term and a potential 
V(9) to the action, leading to the equation of motion 



D0- 



dV 
~d9 



1, 



-R R. 



(7) 



For the action in Eq. © to remain invariant under the 
shift symmetry 9 — > 9 + c, where c is a constant, the 
potential must vanish, rendering 9 a massless field. In 
addition to this, the dynamical theory also adds a stress- 
energy tensor for 9 to the right-hand side of Eq. ([5]) . 



Consider a circularly-polarized, plane GW, propagat- 
ing on a flat Minkowski spacetime, of the form 

Vl(U'') = V R , L cxp {~i y>(t) - kix 1 ] } , (8) 

where t and x % are time and space coordinates, k % is the 
spatial wave vector and /io,r.l is some time-independent 
amplitude. Alternatively, one can work directly with the 
Fourier transform of the waveform 
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where we assume the Fourier modes take the form 
fc R .L,ik(t) = fto,R 1 Le-^W. 

After linearizing the CS modified field equations and 
using Eq. ([5]). the evolution equation for the phase in a 
right/left circularly-polarized basis becomes 



i<j> + 4f 
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where overhead dots stand for time derivatives, k = 
(kik 1 ) 1 ' 2 is the wave number and where A n = +1 and 



A L = — 1. Clearly, the phase associated with the right 
and left polarization obeys different evolution equations. 
Formally, the phase should also have a subscript r,l, but 
we drop it for now to avoid notational clutter. 

Consider now a circularly-polarized, plane GW, prop- 
agating on a Friedmann-Robertson- Walker spacetime, of 
the form 

hnAv, X % ) - Vr,l exp {-i [<j>(rj) - mx*] } , (11) 

where r\ and x % are the conformal time and space coor- 
dinates, and k 1 is the conformal wave vector. We here 
let 4>(rj) be a complex function, such that its imaginary 
part leads to an amplitude time-dependence. Just as be- 
fore, we can alternatively work directly with the Fourier 
transform of the circularly-polarized waveforms 
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where now the Fourier modes take the form h R ^ K (r]) = 

"0,R,L c 

Once more, after linearizing the modified held equa- 
tions in a left/right circularly-polarized basis one finds 
that the phase obeys a similar equation [lj| [21j 
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where % = a' /a is the conformal Hubble parameter, with 
a the dimensionless scale factor that relates conformal 
and physical time (dt = a drj), primes are conformal 
time derivatives and k = (Kin 1 ) 1 ' 2 is the conformal wave 
number. One can easily show that Eq. (fT3|) reduces to 
Eq. ([TU]) in the limit a(rj) — > const. 

In the flat spacetime case, the modified phase can be 
obtained by perturbatively solving Eq. (flUl) : expand the 
phase as 



(14) 
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and 
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with 5(f) <C (f> and assume that 
The GR solution then returns <f> = ■■ ±k. For weak CS 
parity violation (k9 <C 1 and 9 <C 1), the correction to 
the phase becomes 



= iA R , L 7r/ / 9dt = iA R , L 7r.f 



(15) 



where we have defined the GW frequency / = k/(2ir) and 
where the subscript "0" and "s" mean the value at the 
observer and the value at the source respectively. Equa- 
tion P^)) reduces to 



iK^nfoD 6 , 



(16) 



if one freezes 9 = Bq = const., where D is the comoving 
distance to the source. Notice that it is the combination 
D f that determines the magnitude of the modification, 
as already discussed in the Introduction. 



In the cosmological case, the CS modification to the 
phase can be obtained in the same fashion as in the flat 
case. Expand the phase as in Eq. (fT4")). but now <f> and 5(f) 
are functions of conformal time and allowed to be com- 
plex. Use this phase ansatz to series-expand Eq. (|13p . 



assuming that 5(f> <C 
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and 5(f" < (f)'t 



which follows from k9' <C 1. The GR solution then re- 
turns <f)' = ±« — iH.. For weak CS parity violation, the 
correction to the phase becomes 



(17) 



iA R)L 7r / f(t)[9-H9)dt 



iA R-L 7ra /o 



"o 



where we have transformed back to cosmic time using 
9" = a 2 (9 + H9), with the Hubble parameter H = 
a J a = aH. Notice that the GW frequency is formally 
time-dependent due to the Hubble dilution, 2nf(t) = 
k — iH(t), but the second term is subdominant and can 
be neglected; cf. Section IIIII and particularly the Ap- 
pendix [Eq. (|A~T2|) ]. Also notice that if $ = or H = 0, 
we recover the flat spacetime result in Eq. (|16|) . Simplify- 
ing Eq. (fT?]) by Taylor expanding 9 = 9 + 9®t and using 
that a s ~ &o(l — z ) w ^ find 
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where Hq = (d/a)o is the value of the Hubble parameter 
today and we have assumed a matter-dominated Universe 
in an z -C 1 expansion, suitable for GW/GRB event. 
The derivation presented is consistent with the results 
found in [18[ , if one expands the latter in the low-redshift 
limit. This expansion is justified since most LIGO- Virgo 
GW/GRB events are expected to have z<l [42l |. 



C. Amplitude Birefringence 
and Parity Violation in CS gravity 

The imaginary nature of the CS modification to the 
phase leads to amplitude birefringence, i.e., the am- 
plitudes of ft- L and /i R are modified differently leading 
to a suppression/enhancement with propagation. For 
weak parity-violation, 5(f) ^C </>, the exponential suppres- 
sion/enhancement of the right/left circularly-polarized 
GW becomes 



'^'R,L — ^RX e 



ft R ,L (1 



X (19) 



where /i RX is the GR prediction for the right/left 
circularly-polarized GW and <50 R]L is given in Eq. (Til?)) 
or p7|) . Notice that we have now reinstated the subscript 
r.l in the phase. Clearly, the parity-violating modifica- 
tion introduces a change in the amplitude since i i50 RL is 
a real number. One should note here that the correction 



leads to finite CS corrections, as 64> RtIj has all along been 
assumed small. In the Appendix we show that Eq. (TT9")) 
does not signal the presence of an instability and that in- 
deed the CS correction is formally bounded from above 
for all relevant values of the CS parameters. 

The degree of birefringence can be quantified through 
v, introduced in Eq. (fT]). This equation can be rewritten 
as /i RL = ^r,l(1 + A RX i>), which can be inverted for v to 
find 



1 (K K\ i 
For weak CS parity violation, we have 
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where wc have defined <50 FLA t 
gating on flat space and S9 FRV 
propagating on a cosmological background. To measure 
this effect, it is crucial to know what amplitude the GW 
signal would have in GR, which is very difficult with 
ground-based detectors alone, except in the case of bi- 
nary inspirals associated with GRB. 

The analysis presented so far has ignored the pole in 
the modified dispersion relation of Eqs. (flO|) and (fT3|) 
by assuming weak CS parity violation. This assumption 
allows us to expand the pole under k9 -C 1. If this as- 
sumption does not hold, then strong CS parity violation 
leads to a resonance in the GW response function that 
might occur inside the sensitivity band of the GW detec- 
tor. This resonance enhances the constraints one could 
place on gravitational parity violation, as we show sepa- 
rately in Sec. [Viand in the Appendix. 



III. GENERIC PARITY VIOLATION IN GW 
PROPAGATION 

Consider a circularly-polarized, plane GW propagating 
on an FRW background. Let us Fourier decompose this 
wave as in Eq. (fT2"|) . In GR, the linearized field equations 
give the wave equation 



h': 



+ K /Jr,l,k — — 2rt/l aI 



(22) 



The first time derivative of the metric perturbation serves 
as a dissipative term that controls how the wave ampli- 
tude drops adiabatically as the universe expands when ti 
is non- vanishing. 

Parity violation leads to a polarization-dependent 
modification to the amplitude evolution, and thus, the 
parity violating generalization of Eq. (|22|) is the replace- 
ment of 



H 



a 
a 



<5r,l 



(23) 



where 5 RX is a polarization-dependent, generalized scale 
factor. For example, In CS theory, the linearized field 



equations reduce to [see e.g., Eq. (18) in J2l| ] 

\H-\ RX K6"/(2a 2 p 



""R.L.K 



k n. R i 



1 - X^K0'/a 2 
where we see this is identical to Eq. (f22j) 



hi 



(24) 
with the re- 



placement of Eq. (|2"3"|) and S R x = ay/1 — X rl k9' /a 2 . 

Solutions to the wave equation can be obtained by as- 
suming a GW ansatz of the form h HIjjK (rf) = e^ l ^^ v \ 
where 4> is allowed to be complex. In GR, this reduces 
Eq. flU]) to 



'f - k 2 = -2in< 



(25) 



where again we assume 4>" -C {4>') 2 ■ For GW wavelengths 
short compared to the Hubble scale, i.e., k 3> "H, the solu- 
tions are <p' = ±k — iH, to leading order in %/ k separately 
for the real and imaginary parts. We arc led then to the 
plane- wave GR solution h K 



■Xv) 






j±iK,rj 



(26) 



which are the same irrespective of the GW polarity. 

The amplitude-birefringent solutions are obtained sim- 
ilarly through Eq. (|23[) . The complex phase then becomes 



±K?y- 



i S' 



R,L j / 

'Jrx 



±k+i In 






(27) 



This results in a plane- wave solution of the form [2l| 
h^, K (v) = ^4e^ 



SkAv) 



(28) 



This solution is subject to the assumption k 3> <S R L /S RX , 
which is satisfied in CS gravity provided 6' varies on 
scales much longer than a gravitational wavelength, 
i.e., 9" /9' <C K. Of course, this assumption does not hold 
in the immediate neighborhood of the singular point of 
the wave equation at <S R>L = 0. This singular point is the 
source of the pole found in the previous section and it is 
addressed separately in Sec. IVland in the Appendix. 

We can now explicitly write the propagator in the form 
of Eq. (UJ. GW polarizations at conformal time r\ are 
given in terms of r\ s near the source as in Eq. (fTJ) with 
the parity violating measure 



v 



1 a(v) 
2a(r) s ) 



Sr(Vs) Sl{Vs) 



Sr(v) Sl(v) 



(29) 



Wc note that v reduces exactly to the defini- 
tion introduced in Eq. ((2H)) when one sets <S R>L = 
ay/1 — Ar L k6" / 'a 2 . 



IV. A GW/GRB TEST OF WEAK PARITY 
VIOLATION 



Let us now derive the effect of amplitude birefringence 
on the GW SNR of a GW/GRB event, as summarized 



by [3J] and references therein. The single-detector GW 
power SNR p after matched filtering is given by the inner- 
product 



(h\h) = 4 / df 



Mf)\ 2 
s h (f) ' 



(30) 



where Sh is the one-sided power spectral strain noise 
density of the detector. The quantity h(f) is the 
Fourier transform of the GW strain measured by a de- 
tector, i.e., the response function h(t), which in a parity- 
preserving theory like GR is given by 



h(t) = F+h+(t)+F x h x (t), 



(31) 



where /i+, x are the GR plus or cross-polarized waveform, 
while F+ and F x are beam pattern functions that depend 
on the sky position and a polarization angle ip (which 
does not concern us here since we deal with nearly circu- 
lar polarization). We use overhead hats to remind our- 
selves that these expressions are valid in GR only. In 
the Fourier domain, the GW power of a GR signal is 
approximately 



IM/)P 



Fl\h + \ 2 



Ft\h, 



(32) 



where h+ tX are the Fourier transform of the GR plus or 
cross-polarized waveform. We have neglected the F + F x 
term here, since it vanishes upon integration in the SNR 
calculation. 

For GW signals from binary NS inspirals observed by 
ground-based detectors, GR predicts that the frequency- 
domain signal is given to a good approximation by 



Hf) = Af 



-7/6 e i*(/) 



(33) 



where the phase *£(/) does not concern us here (except 
indirectly as it yields precise mass measurements which 
help determine the distance). We are interested in the 
amplitude coefficient 



A 



~M 5/6 



96tt 4 / 3 D 



(34) 



where M. is the "chirp mass" and D is the distance from 
source to detector. We neglect corrections to the wave- 
form due to propagation on a cosmological background, 
as most GW/GRB events will be at low redshift. One 
can encapsulate the sky position, polarization, and incli- 
nation dependence in the factor 



2 c i 



Q 2 = Fl(l + C 2 ) + 4F 2 £ 



(35) 



where £ is the cosine of the inclination angle between 
the orbital angular momentum and the line of sight from 
source to detector. 

A GR binary without orbital precession has a squared 
SNR of 



M 5 / 3 



247r 4 /3 D 2 



Q 2 ir, 



(36) 



and thus an SNR of 



1.3 



M 



1.2 M Q 

h 



4.7 x 10 44 s 



5/6 



1/2 



300 Mpc 
D 

F 2 + + Fl 
2/5 



1/2 



(37) 



where in this last relation we have averaged over angles 
and assumed an inclination angle such that £ 2 = 1. The 
noise integrals are defined as 



I P = 



df 



f 



-p/3 



Sh(f) 



(38) 



1, 



Note that Q, and therefore p, is maximized for £ 2 
i.e., looking down the angular momentum axis. 

In the limit of strong signals and weak Gaussian noise, 
the standard deviation of a single-detector GW ampli- 
tude measurement is 



1A 
A 



1 
9 



(39) 



With a GW measurement alone it is difficult to convert 
this into a distance measurement because sky position, £ 
and ip are not easily measured for brief signals such as 
inspirals. A GW/GRB coincidence provides the needed 
information. GRBs with observed afterglows have negli- 
gible sky position errors. If the GRB jet is tightly col- 
limated, £ 2 = 1 and ip is irrelevant. In that case, the 
best-case GW distance measurement error 



D 



1 
9 



(40) 



is achievable [34j. Recently, [35| pointed out that the 
above might be an overestimate of the accuracy due to 
degeneracies between D and £ 2 . However these degen- 
eracies can be lifted by observing numerous GW/GRB 
events with a network of non-aligned detectors such as 
LIGO- Virgo. 

The single-detector GW SNR is altered in parity- 
violating gravity as follows: Combining Eqs. (fTJ) and (|31[) . 
the strain observed by a detector is seen to be 

h = F+(h+ + ivhx) + F x (h x -~iv~h+). (41) 

and the GW power in the Fourier domain is 

\h\ 2 = \h\ 2 + 2 v (F 2 + F 2 ) lm(h x h* + ) + 0(v 2 ), (42) 

where again we have neglected the cross terms of the 
beam pattern functions, as these vanish upon integration. 
This relation then changes p 2 from the GR value p 2 by 
the replacement 

Q 2 -> Q 2 + 4v (F 2 + F 2 ) £ (1 + e) + 0(v 2 ), (43) 
in Eq. (J3BJ, with Q 2 given by Eq. §E§. 



Assuming that v -c 1, we find that the change in rela- 
tive SNR due to weak parity violation can be expressed 
in terms of geometrical factors times a weighted average 
of the model-dependent factor v: 



P 



= 1 



(Fi+Fi)z(i+e) 



I r ^/~ 7/3 

I7J0 f s h (f)- 



(44) 
(45) 



For a perfectly collimated GRB (£ 2 = 1) this reduces to 
the sky-position-indcpcndcnt first-order result 



P 



1 + A RJ 



(46) 



where the factor of A R L arises because £ = ±1 depending 
on the orientation of the face-on source and we have asso- 
ciated a source with £ = +1 with a right-polarized wave. 
For a weak CS modification to gravity in an FRW back- 
ground, v reduces to Eq. (f2Tj) and the SNR correction 
reduces to 



- = 1 + A RiL 7T S9 FRW — . 
P h 



As in Eq. (f!5]l. in the limit of slowly- varying 9, Eq. (|4"7|) 
reduces to 



4 = i-;k,*Wflo-|r]£, ( 48 ) 

P \ H / h 



where we have used that D ps z/i7 f° r zd- 

In the best case, p is a Gaussian with variance unity, 
and thus a (lcr)-bound on generic parity violation trans- 
lates into 



(v)<l/p. 



(49) 



For weak CS parity violation with 9 = 9o, this test would 
allow one to place the constraint 



i<bi<i£i. 

IT I4 Z p 



(50) 



With Eq. (|36p . Eq. (J5U| can be rewritten in terms of 
observables as 



K'\ > y/l^'y/n + nHoMV*-^. 



(51) 



For a double NS GW/GRB, averaging over sky 
location (such that Ff = F 2 = 1/5), H a = 
70.5 km s _1 Mpc -1 , and a single interferometer with 
the present estimate of the "advanced LIGO" fully- 
commissioned noise performance [43], [44j , the threshold 
is 



|0 O | < 1200 km 



1.2 M< 



5/6 



2/5 



F| 



1/2 



(52) 



Therefore, the 95% confidence threshold for parity viola- 
tion (compatible with a 2cr constraint) translates into the 
bound \9o\ < 2400 km. Note that the dependences on z 
and p have cancelled out. However, p does matter in that 
z-independent factors which improve p, also improve the 
threshold on 9q, up to a maximum p set by systematic 
errors (see Sec. lVip . 



V. RESONANT MAGNIFICATION 

The analysis presented above neglects a resonance in 
the parity violating interaction [see e.g., Eqs. (fi"0|) . (fT3|) 
or Eq. (|29p ]. As is apparent from the ubiquitous denom- 
inator 1/<S R L in v, certain wavelengths will be strongly 
amplified. Such an amplification leads to a finite increase 
in the amplitude correction, as it is shown in the Ap- 
pendix. Without specifying a functional form for <S R-L , 
it is nearly impossible to study such resonant behavior, 
which is why in the rest of this section we assume S RL is 
given by the CS prediction in Sec. IIIII 

In CS gravity, the resonant frequency f p and resonant 
conformal time rj p are such that <S RL = 0, which reduces 
(47) to the condition n p = A RlL Op/#o or i n terms of cosmic 



time 



fp,o 



Ar,l 

2tt9 



(53) 



Nearly the same frequency band will be resonant for all 
GW sources, because the resonance condition depends 
only on 9$ at the detector. 

The estimates of Sec. IIIII for amplification are ade- 
quate, provided the GW frequency is sufficiently far from 
this present-day resonance / P) o! the resonant frequency 
at the source f p>s = A RL /27r6> s ; and any frequency in be- 
tween, as these frequencies have generally been resonant 
in the past. As described in the Appendix, the strongest 
amplification will occur very close to the present-day res- 
onant frequency / Pi o, with a bandwidth of 5f p w (fp) 1 , 
which reduces to Sf p = [9/(2n9 2 )] 1 / 2 . Assuming that 9 
evolves on cosmological scales 9 = 0(H9) implies that 
5f p = 0[{H/9) 1 / 2 ] = O(10- s Hz). 

Similarly, the time period during which the signal is 
resonant can be computed from the inverse of the reso- 
nance bandwidth, or equivalcntly from the rate of change 
of the resonant conformal time r\ p with respect to con- 
formal wavenumber k: 1Z = drj p /dK p . Let us define the 
squared amplitude J^l = 5r L j which is a function of 
two independent variables: r\ p and k p . Requiring that 
dJ- BL = leads to 



drj p 



Uj~ R l/ UK p 

dT K ±/dr) p 



and rewriting this equation slightly we find 



n 



d 
dn 



T' 



(54) 



(55) 



where we have used the fact that .F R]L = at resonance. 
Although the resonance bandwidth is small, if a reso- 
nance is presence in the sensitive band of the GW detec- 
tor it will greatly enhance the constraints of the previous 
section. In what follows, we compute the enhancement of 
the simple GW/GRB test discussed in Sec.|lV]if the res- 
onance is in-band. We begin by assuming that the GW 
signal detected is parity violating and one searches over 
the data with parity-violating templates. Such consid- 
erations lead to an optimal constraint related to the in- 
trinsic SNR of the parity- violating templates. We follow 
this analysis by considering a parity- violating signal that 
is filtered with a GR template. By doing so, one incurs 
a fundamental error in the determination of parameters 
induced by the assumption that the data contains a GR 
template (see e.g., |45( for a more detailed discussion of 
fundamental bias in GW astrophysics). 



A. Optimal Strategy 

1. Signal-to-Noise Ratio and an Optimal Test 

The presence of a resonance in the frequency band 
where searches are performed leads to an amplification 
of parity violation. Such an amplification can lead to a 
stronger constraint than the one we considered in Sec. lIVI 
Of course, the constraint is largest when the search is 
performed with templates that account for such parity 
violation. We consider such cases in this subsection. 

Let us then assume that some of the detected frequen- 
cies are resonant and that one correctly identifies all sig- 
nal power associated with the event to the correct parity- 
violating template. The SNR then becomes 1 
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(57) 



where we have defined F R , L = (F + + A R , L F x )/\/2 and 
X R>L = (j r R , L , s /.7 r Rl i,,o) 1 i an d where we have assumed 
h + and h x are uncorrelated, so that their integrated 
product vanishes. Recall that hatted quantities cor- 
respond to the GR expectation, so that h +x are the 
Fourier transform of the GR, plus or cross-polarized 



waveform. Notice that this expression reduces to the 
GR result when X B = X^ = 1 and to the expression in 
Eq. (03J in the linear approximation X RjL ~ I + A RL w. 

The above expression can be simplified somewhat if we 
assume the source is at £ = ±1. We then find 



7^5/3 



127r 4 /3 D 2 



(f 2 +f 2 ] 



x r f r V3 [(Xl + Xl) ± (Xl Xt)] .(58) 

JO °h 

The beam-pattern functions then only modify the overall 
coefficient of the SNR, and thus we can rewrite the SNR 
in the more familiar form 

p - = 2A 2 ^ ^r 7/3 [{Xl + Xl) ± (Xl Xl)] , 

JO &h 

(59) 
Given a single detector, we cannot distinguish between a 
purely right, or purely left-polarized signal, so effectively 
we measure either p R or p L , where we have defined 



pi., 



A 2 



1 



t-7/3 



J~ R.LWSI k) 



and the amplification factor is 

1 r°° f-r/3 

Qn,, 2 = T / dfl—- (1 + z) 
h Jo i/il/J 
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J^rAvo^) 



(60) 



(61) 



We have here used that p = 4 A 2 Ij/(\. + z) as shown in 
Eq. (f36|) . but we have included a factor of (l+z) -1 to map 
between comoving and luminosity distances. This result 
generalizes the linear, weak-amplification expression of 
Eq. dHJ. 

Let us now show that the above reduces to the results 
of Sec. IIVI in the weak parity violation limit. Assuming 
the resonant frequencies are very large (i.e., ko9 ^C 1) and 
that the source is at low redshift DHq <C 1 , the resonant 
factor becomes to leading order 



J" RX (?/ S ,fc) 



FnAvo,k) 



1 - 27rfDX R , L e H - e 



(62) 



We recognize the right-hand side as simply l + 2A RL i> and 
inserting this back into the SNR we find 

i r°° f-7/3 
P 2 - P 2 T dfJ—-(l + 2\ R ^v), 

h Jo bh(j) 



1 

h Jo 
f? (1 + 2A R , 



df 



f-r/3 
SUf) 



1 r°° f-r/3 



(63) 



1 Technically, the one-sided SNR integral arises from a more gen- 
eral, two-sided overlap integral l46l . acting on a real- valued 
(i.e., single polarization) strain signal h+(t). Even though the 
J-r,l is not symmetric in /, the Fourier transform can be folded 
together using that J i r(/) = J-"l(— /)■ 



where (v) was defined in Eq. (|45|) . Notice that this ex- 
pression reduces to Eq. (|46[) , upon linearizing the square- 
root. 

A test of parity violation then presents itself. Let us as- 
sume that a GW detection has been made coincident with 



10 



a GRB event at a known distance, so that the chirp mass 
and p (the GR expectation) are unambiguously known. 
Then, parity violation can be distinguished from a (2a)- 
statistical fluctuation if |p R>L /p — 1| > 2/p [cf. Eq. (|50"|) ]. 
If such a violation is not observed, one can rule out points 
in the (9,9) space. In particular, (9q,9q) values that lead 
to 



|S R , L -1|>-, 
P 



(64) 



can be excluded; inversely, (9o, 9q) values that lead to the 
opposite relation cannot be distinguished from pure GR 
with a single observation. 



2. Numerical Implementation 

The implementation of the test in Eq. (|64p requires 
the calculation of the amplitude enhancement factor <7 R , L , 
which in turn requires knowledge of the ratio of J RL at 
the source and at the observer location. The latter can 
be rewritten as 



(1 + z) 



■TR^lls] fc) 



J^rAvo^) 



1 - f/(a s f P ,s) 



1 - f/(aof P ,o) 



(65) 



where ao/a s = 1 + z and we have defined the resonant 
frequencies at the source f PiS and at the observer / p .o- 
The resonant or pole frequency at the observer is simply 
fp,o = ^r,l/(27t0 P! o): while that at the source can be 
approximated as 



dt 



fp,8 — fp,0 + -77<>t — fpfi - 



D 

2wK 



f p ,oH Q D, (66) 



where we have used that / = k/(2tto), f p = n' /(2-ko 2 ) 
KpH/(2Tra 2 ), k' = Hr 1 and St = —D. Combining f p 
with the appropriate redshift factors that arise in a s : 
(1 — z)ao, we find 



&s Jp,s ~ Jp,0 
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2-kTI 



OiKfiHf 



(})• 



(67) 



where the last term of Eq. ([67]) has cancelled. 

The calculation of the enhancement factor then re- 
quires knowledge of 1Z, which following Eq. (|55[) reduces 
to 



K 





-I -2 



(68) 



Using further the fact that k p ,o — f p .o/(2ir) 
A R x /(47r 2 0o), we can simplify the above expression to 



n = 



i 



2irf pfi H Q (l-q) 
9q/(9qH ) , 



(69) 

(70) 



where the dimcnsionless parameter q is 0(1) for slowly- 
varying 9. Using this expression and in the limit z ~ 
H n D <C 1, Eq. ([67]) reduces to the extremely simple re- 
lation 



f P ,s « f p , [l + z(l -q)] 



(71) 



In the local universe or for slowly-varying 9, the 
present-day and source-frame resonant frequencies are 
nearly identical. More precisely, the first-order approx- 
imation used in Eq. (|67p to relate these frequencies im- 
plicitly requires z\(l — q)\ ^C 1, or 



zti 



H a 9 Q \ « \H 6, 
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(72) 



corresponding to the requirement that 59 FKW be small. 
The analysis presented here does not formally apply when 
the pole frequency is rapidly varying, i.e., when Eq. (|T2[) 
is violated. However, even then, non-resonant and reso- 
nant amplification will gcnerically occur and one should 
still be able to constrain gravitational parity violation. 
The results we present below, however, are not valid in 
this region. 

The integrals in Eq. (f50|) and (|rJT|) are evaluated with 
the estimated NS-NS optimized spectral noise density 
of [43| and numerically converting them into Riemann 
sums and the Fourier integrals into Fourier series. In 
doing so, we must choose a frequency discretization df, 
which is related to the observation time T via df ~ 
1/T. Full coherent recovery of the amplified waveform 
clearly corresponds to an unrealistically long coherently- 
integrated data stream: for a T = 1 year long integration, 
df = C(l/yr) ~ 10~ 8 Hz. In practice, due to detector 
dropouts and nonstationarity of the noise, the longest 
integration times are 2048 s |47j . imposing a minimum 
frequency discretization of df ~ 1/T ~ 5 x 10~ 4 Hz. 
A search strategy based on these limited data intervals 
will always recover a somewhat smaller resonant contri- 
bution to the total SNR than if one where to integrate 
for larger T, which then implies less limiting constraints 
on (9q,9o), assuming the same SNR threshold. However 
the dependence on df is generally quite weak. We have 
verified these statements by numerically computing the 
integrals in Eq. ([60]) and (|BIJ both with df = 10~ 8 Hz 
and df = 1/2048 Hz and checking that the changes to 
Fig. [2] (described below) are hardly visible. In what fol- 
lows, we shall adopt a realistic frequency resolution of 
df = 1/2048 Hz. 

The integration of the resonant amplitudes arc also for- 
mally ill-defined when the integrand sweeps through the 
resonance / — ¥ / p ,0i as then ^r.l^Oi^) vanishes. Such 
divergent behavior can be cured by regularizing all in- 
tegrals, i.e., by cutting out a region of size e about the 
resonant frequency at present (f P ,o)- Such a regulariza- 
tion is valid because the SNR depends logarithmically on 
the regulator, as we show explicitly in the Appendix. We 
choose here e = Sf p , but we have checked that choos- 
ing e = df = 10 -4 Hz, for example, docs not visually 
affect any of the plots. This regularization is allowed 
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FIG. 2: The blue shaded region corresponds to values of 
(9o,6o) for which parity violation is too small to be mea- 
sured with a coincident GRB/GW event. The region outside 
this blue-shaded area can be measured or excluded given a 
single, GW/GRB event at p — 10. The boundaries of the 
blue-shaded region correspond to the best (#o, #o) constraints 
one could place. The thick solid line corresponds to the val- 
ues of (8o,8o) for which there is no parity violation (Q n = 1), 
and thus it represents an absolute theoretical limit to the con- 
straints one could place on CS amplitude birefringence. The 
dotted lines bracket the values of (0o,0o) that are allowed 
by Solar System constraints. Consequently, the regions out- 
side the dotted lines (not shaded in the figure) have already 
been excluded by Solar System tests. The green-shaded re- 
gion shows where the linear, weak-parity approximation of 
Sec. IIVI is valid, while the gray-shaded region shows where 
the linear cosmological expansion of Eq. (|67p breaks down. 
Our calculations suggest that LIGO/Virgo will be able to see 
or rule out the entire area outside the blue-shaded region, ex- 
cluding a much larger region of parameter space than Solar 
System constraints. 



within the context of linear propagation theory because 
the resonant amplification enhances the waveform only 
to hf/Sfp <C 1. Even though the enhancement factor is 
formally large, f/Sf p ~ O(10 10 ), the overall correction is 
clearly small as typical GW strains are h = 0(1O~ 22 ). In 
the Appendix, we show in more detail that regularizing 
these integrals induces a negligible error in the SNR. 



3. Results and Discussion 

Figure [2] shows numerical results for the implementa- 
tion of the test devised in Eq. (|M|) . The blue-shaded 
region, labeled GR consistent, corresponds to values of 
(0o, 9q) for which we cannot distinguish between GR and 



CS gravity, i.e., values for which \Q R — 1| < 2/p = 2/10. 
This region is generated assuming D = 600 Mpc and 
p = 10 as a fiducial threshold, which could occur in ad- 
vanced LIGO given a NS-NS merger at an optimal sky 
location. Given a coincident GRB/GW detection that 
is non-parity violating, LIGO/Virgo will be able to ob- 
serve or rule out the area outside the blue-shaded region 
in the figure, which is why we have labeled it measur- 
able/excludable. Observe that the blue-shaded region in 
this figure corresponds to the non-shaded region in Fig.Q] 

The boundary of the blue-shaded region in Fig. [2] de- 
notes the best bound that one could place on (0o,0o) 
given a coincident GRB/GW detection consistent with 
general relativity. This boundary corresponds to values 
of (0 O , O ) for which \Q R - 1| = 2/p = 2/10. The Solar 
System constraint 6q < 2000 km with 0o = is out- 
side the plotting range of this figure. One can generalize 
or improve on the Solar System constraint via the map 
|$o| —I \Qo — &o/Ho\, assuming the entire right-hand side 
can be treated as a constant, evolving on timescales much 
larger than the orbital velocity of the LAGEOS satellites. 
Doing so, the improved Solar system constraint becomes 
|#o — 0o/#o| < 2000 km. We have plotted this constraint 
with two diagonal dotted lines in Fig. [2] The region be- 
tween these dotted lines is consistent with GR and can- 
not be ruled out, while the area outside of the dotted 
lines can be excluded with the improved Solar System 
constraint. 

We can now compare this generalized Solar System 
constraint (the dotted lines) to the possible GW/GRB 
constraints (the boundary of the blue-shaded region). 
Figure shows that the latter largely is inside the for- 
mer when 0o > 0, which implies GW/GRB constraints 
would be better than the generalized Solar System con- 
straints. When 0o < the GW/GRB constraint is not 
as good with a single detector, as CS parity violation 
would then be non-resonant in Q n . Such an issue is 
circumvented if either (a) two GW detectors are used 
to extract both h R and /i L for a non-optimally aligned 
source |£| ^ 1 (we discuss this in the next paragraph) 
or (b) two sources, one right- and one left-handed, have 
been observed. Furthermore, if we assume 0o = as in 
the standard Solar System constraint, or more generally 
whenever |0o| < 100 Ho km s _1 , our results suggest that 
a GW/GRB detection could constrain O < 100 km. This 
bound is 20 times stronger than the current Solar Sys- 
tem constraint (0o < 2000 km). Such enhancements in 
the GW constraints relative to those presented in Sec. IIVI 
are clearly due to the CS resonance in the SNR. 

Resonance amplification for a right-circularly polarized 
GW only occurs when the integrand of the amplitude 
factor Q R contains a pole, which leads to the logarith- 
mic divergence discussed in the Appendix. For 0o < 0, 
the amplification of right-handed signals is suppressed as 
the pole moves to negative frequencies; constraints from 
right-handed measurements on O , O arc therefore not as 
strong. The situation reverses, however, if the GW is left- 
polarized, as then the C/ L regions are identical to the Q R 
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ones in Fig.[2]but with (6 o ,0q) -» (— 0o, — 0o)- Therefore, 
the two-detector detection of a a right and a left-polarized 
GW would allow us to constrain (60,60) equally well for 
both positive and negative 0o- More precisely, with such 
a detection one would be able to constrain |0o| < 100 km 
for any |0 O | < 100 H km s" 1 . 

An interesting feature of CS parity violation not ap- 
preciated with Solar System constraints is that for cer- 
tain values of (0o,0o) amplitude birefringence can actu- 
ally vanish. In some sense, for these values of the deriva- 
tives of the CS coupling, the CS correction conspires with 
the cosmological expansion to lead to no parity violation. 
Such values of (6o,6o) then constitute the ultimate the- 
oretical limit to which CS gravity can be constrained. 
These values of the derivatives of the CS coupling are 
such that Gr,l = 1, which implies p = p. This degener- 
acy occurs when / Pj o = fp.s or simply when #o = Hq9q, 
which is simply a straight line in the (0o, 00 ) phase space. 
We have plotted this linear relationship as a solid diago- 
nal line in Fig. [2] Observe that this line is always inside 
the blue-shaded region, as it should, as the latter rep- 
resents precisely the region where a GW/GRB measure- 
ment does not have the sufficient accuracy to see parity 
violation given an SNR of 10. 

Although resonant amplification generally occurs for 
any prescribed form for 9, the Taylor expansion used in 
Eq. (|66| to relate source quantities to observer quantities 
breaks down for certain values of (9q, 9q). More precisely, 
when (0q, 0o ) is such that |^(1— q)\ = 0(1) or equivalently 
\H 9o — 0o | = 0(H 9o/z), the pole frequency changes 
rapidly and the Taylor expansion of Eq. (|66p is not valid. 
One might then worry that our analysis breaks down in 
this region and the bounds quoted earlier in terms of the 
boundary of the blue-shaded region would not be valid, 
as higher-order terms in Eq. (|66p would be needed. To 
get a sense of the (0o, 0o) values for which this happens, in 
Fig.[2]we have shaded in gray the region where \z(l—q) \ > 
0.5 and labeled it rapid pole variation. Observe that this 
region does not overlap the blue-shaded region for 9 > 0, 
which explicitly shows that the bounds proposed above 
still hold. Of course, one can shrink this gray region by 
Taylor expanding Eq. (|66| to second-order if needed, at 
the expense of introducing a new parameter 0o- 

The numerical results we have presented here do re- 
cover the linear approximation [i.e., Eqs. (|48jl and (|18j) ] 
for the (0o,0o) values where this approximation is valid, 
i.e., where there is no resonance in the detection band. 
The source and observer resonances are out of band pro- 
vided either f Pt o < — (/) and f Pt o > 300Hz (for the sup- 
pressed frequency being negative or out of the most sen- 
sitive band) or f PtS > 500Hz and f p . s < — (f) (when the 
resonance is well out of band or negative). In order to 
visualize this in the context of the bounds placed by the 
boundary of the blue-shaded region in Fig. [2 we have 
shaded in green the linear approximation region and la- 
beled it linear CS . Observe that the linear approximation 
breaks down for all physically interesting and accessible 



0o > 0: a resonance generically appears inside the LIGO 
band. Notice that this is not the case for LISA, as this 
instrument is sensitive to much lower frequency GWs, for 
which the resonance is always out of band if one saturates 
0o with the Solar system constraint. 

Finally, we have checked that the constraints presented 
here are formally independent of distance to the source 
for (0o,0o) in the linear regime, which is consistent with 
the observations of Sec. lIVI This occurs because although 
CS amplitude birefringence is proportional to distance 
[see e.g., Eq. (HE])], the dominant source of error also 
grows with distance. Therefore, to leading order, any 
single, short GRB observation consistent with pure GR 
translates into upper limits on 0o — 9 /H$ that are in- 
dependent of distance and depend only on the absolute 
detector sensitivity along the line of sight [Eq. (f5"2")) ] 2 . 
On the other hand, in the non-linear regime, where reso- 
nant behavior is present and the linear approximation is 
not valid, the blue shaded region in Fig. [5] does depend 
strongly on source distance. 



B. Suboptimal Strategy 

The previous analysis assumed the entire signal, in- 
cluding its narrow resonantly amplified contributions, 
arc correctly identified through an appropriate, parity- 
violating template and summed in phase to extract the 
relevant signal power p. In other words, in the previ- 
ous subsection we assumed a follow-up of each detected 
signal with a family of CS parity- violating, resonantly- 
amplified templates possessing two additional parameters 
(/p,sj fp,o) or equivalently (0o, 0o). In this section, we in- 
vestigate constraints on parity violation if all matched 
filtering is performed with parity-preserving, GR tem- 
plates. 

Given a GW detection, matched filtering with GR tem- 
plates will return a set of best-fit parameters A GR that 
maximize the SNR p max . In the absence of noise, the 
values of (0o, 0o) that lead to p = p max arc simply those 
on the diagonal of Fig.[2J i.e., 0o = Hq9q. In the presence 
of noise, however, random fluctuations in the amplitude 
can resemble a resonant amplification. Therefore, one 
can search for the largest values of (0o,0o) that lead to 
a p that is consistent with p max . Such values of (0o,0o) 
then serve as an upper limit or constraint on the magni- 
tude of parity violation in the signal. 

Such a study leaves out the fact that a GR matched 
filtering calculation unavoidably induces a fundamental 
error in the estimation of GR parameters, if the signal is 



2 As in the linear case, the expression (v) < 1/p mislcadingly sug- 
gests distance-dependent results. Areas in Fig. J2I only apply to 
self-consistent choices for p and D. For example, the contour for 
\Q — 1| < 2/10 at D = 600Mpc corresponds to either twice as 
sensitive an instrument or 2 6 ' 5 times larger a binary chirp mass. 
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parityviolating. This is a manifestation of fundamental 
bias [45| . In other words, the true physical parameters 
of the signal A = A GR + SX, where SX is the aforemen- 
tioned error. We shall here ignore this error and leave an 
investigation of fundamental bias to a future publication. 
As in the previous subsection, given a source with |£| = 
1 and a single detector, we cannot distinguish between a 
purely right or purely left circularly-polarized signal, and 
thus, the SNR p we would detect is 
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Notice that we have here crudely assumed the resonantly 
amplified waveform is projected onto a pure GR signal 
model with identical parameters. To improve the realism 
of this calculation, one could maximize over time and 
phase of coalescence, but we shall not do this here. 

Unlike the intrinsic SNR of the previous subsection, 
the integral defining p requires no regularization. This is 
because G B L differs from C/ R L in that the former contains 
a square root, while the latter contains an absolute value. 
This square root leads to a well-defined integral, as we 
explore in more detail in the Appendix. All other details 
of the numerical implementation are the same as in the 
previous subsection. 

When resonant frequencies are very large (i.e., in the 
linear regime, where constraints on 8' are weak), then 
Eq. ([71]) reduces to 
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which is precisely half the amplification seen using an 
optimal signal model in the previous subsection. We have 
computed figures similar to Fig. [5] and we indeed find 
that the constraints quoted in the previous subsection 
deteriorate at most by a factor of 2. Of course, when 
one maximizes over extrinsic and intrinsic parameters, 
the match will increase, and thus the constraints will 
improve. 

VI. EXTENSION AND DISCUSSION 

In this section we discuss estimates of various factors 
that improve the order of magnitude results from the 
previous section. We also discuss limiting noise sources 
that could potentially deteriorate the constraints. 

A. Constraints from single short GRBs 

A network of Nd GW detectors using coherent data 
analysis is more sensitive than one detector by a factor 



y/Nd in p, neglecting differences in orientation and noise 
spectrum. The four-detector LIGO- Virgo network thus 
gains about a factor of 2 for a typical sky position, though 
the full beam pattern is complicated |4(}|4]|. Moreover, 
this network also is more robust than a single detector if 
the GRB is viewed off-axis [HHi]. 

A BH/NS merger takes advantage of the mass de- 
pendence of the threshold (fSTj) . Ai 5 ^ e , which goes as 
the square root of the total mass. Thus a 30 M Q BH 
(M. = 4.7 Mq) would improve the threshold further by a 
factor of about 4. Above 30 M Q the improvement drops 
because the system starts merging before the frequency 
has swept through the LIGO band and because our anal- 
ysis assumes no spin-orbit precession. If the BH is rapidly 
spinning and that spin is misaligned with the orbit, pre- 
cession will lead to a time-dependent left/right polariza- 
tion ratio, reducing p from the on-axis case and com- 
plicating the analysis. Spin-orbit misalignment could be 
large for dynamically-formed binaries. In BH-NS binaries 
formed from isolated binary evolution, spin and orbital 
angular momenta should usually be aligned and preces- 
sion should be relatively unimportant; see e.g., Fig. 6 in 

Observations and population arguments suggest that 
GRB associated with NS/NS mergers have a narrow but 
nonzero range of opening angl es; p opulation arguments 
suggest narrow jets 0^2 = 0.02 [5(|, while the lack of jet 
breaks in af terg lows suggests that some jets may be as 
wide as 25° pCf (which is still £ w 0.9). The finite colli- 
mation angle of the GRB will lead to a distribution of £ 2 
which will degrade the threshold if a ^ ~ 1/ p since a^z 
adds roughly in quadrature with ap/D. Additionally, 
because £ 2 < 1, these unknown angles can only reduce 
p, not increase it, introducing a weak negative bias to- 
wards smaller amplitudes. Without independent direct 
or a priori constraints on £, any short GRB amplitude 
measurement by itself can only constrain the particular 
combination A R , L (v) + (A RX £ — 1), as 
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where we have explicitly expanded Eq. (|44[) for a nearly 
face on, right- or left-handed source £ w A R L and used 
that />(£) ~ /3(Ar, l )[1 + A R>L (£ - A R>L )]. In practice, how- 
ever, the LIGO- Virgo network (or any network of dif- 
ferently aligned detectors) is sensitive to both GW po- 
larizations, and by projecting directly onto a circularly- 
polarized basis, this network can directly constrain 1 — 
£ ~ 0(1/ p). [A similar effect was observed in Seto [33[ 
and Dalai et al. [34| , where the LIGO- Virgo network (or 
any network of differently aligned detectors) can reduce 
the systematic errors in distance measurement to short 
GRBs induced by wide jets.] The net effect of a nar- 
row distribution in £ is therefore simply to weaken any 
single-observation constraint on (v) by roughly a factor 
of 2, from 1/p to 0(2)/ p. 
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Another source of error is due to mismodcling of 
the signal because of the neglect of higher harmonics 
present for example in PN amplitude corrections j52| . 
Such corrections appear in the stationary-phase of the 
Fourier transform as modifications to the amplitude of 
the form A ->■ A (l + Ci/2" 1/3 + (m 1/3 + C3/2"), where 
u = 2nMf is the reduced GW frequency and Cn are nth 
order PN amplitude coefficients, which depend on the 
mass ratio, the inclination angle and the polarization an- 
gles. Notice that the 1.5PN correction, £3, is proportional 
to /, which immediately suggests a strong-correlation 
and degeneracy with CS amplitude birefringence. This 
degeneracy should be affect the constraints quoted above 
much because the mass ratio is primarily determined by 
the phase at 1PN order, while < v > does not enter the 
phase at all. A full study of such possible degenerate 
behavior would require a full Fisher or Bayesian analysis 
that we shall carry out elsewhere. 

The electromagnetic determination of the distance to 
the GW/GRB also has errors. The latest augmented 
WMAP results [&| imply H = 70.4±1.3 km s" 1 Mpc" 1 
at 95% confidence. The intrinsic uncertainty in the Hub- 
ble constant translates to a relative distance uncertainty 
of roughly SH = 1.3/70.4 < 1.8%, far below the intrinsic 
SNR uncertainty 2/p for all but extremely bright sources 
(p > 108). Peculiar velocities of host galaxies, internal 
orbital velocities within galaxies and galaxy clusters, and 
kick velocities of binaries due to supernova kicks should 
all be well under 1500km/s ~ 0.005c (54J, correspond- 
ing to a distance error of about 6% in a Hubble redshift 
of z = 0.08. Only for extremely bright sources or very 
close double NS binaries will non-cosmological velocities 
produce a noticeable error. 

A more significant and less well controlled distance er- 
ror arises through the redshift determination process it- 
self. Currently short GRB distances can be estimated 
only through the redshift to electromagnetically asso- 
ciated hosts, and therefore depends critically on cor- 
rect host identification |55l457| |. For moderate to high 
redshift sources, particularly where short GRBs are as- 
sumed ejected from their putative host galaxy, associa- 
tions and therefore redshift measurements could be sig- 
nificantly confused J58(. Fortunately, for the distances 
to which even advanced gravitational detectors are sen- 
sitive, relatively few galaxies can conceivably be asso- 
ciated with a coincident gravitational and electromag- 
netic signature (i.e., adopting a fiducial galaxy density 
n gai = 0.01 Mpc -3 , the number of galaxies in the elec- 
tromagnetic sky-position angular extent 9 and radial ex- 
tent D error cone is of order n ga iD 3 ird 2 /3 = 0(1) for 
D = 600Mpc and 8 ~ larcmin). We therefore assume 
any future coincident short GRB and GW signal arising 
within z < 0.4 will have a correct host association and 
therefore a reliable redshift measurement. 

Let us then summarize the enhancement and deterio- 
ration factors we have just discussed. A network of Nj 
detector enhances the bound by a factor of y/Nn, while 
a BH/NS event enhances the constraint by a factor of 



4. On the other hand, an unknown jet width deterio- 
rates the constraint by a factor of 2, while a suboptimal 
search strategy also worsens the bound by a factor of 2. 
Putting all of these together, we expect a modification of 
the constraints quoted earlier by a factor of l/y/Nj, for 
a NS/BH event or A/y/N2 for a NS/NS event (i.e., from 
(v) < 1//5 to l/(p^/Nd) for a NS/BH event, where p is 
the single-detector SNR). 

Putting all of this together, we arrive at the follow- 
ing constraints. A single detector with advanced LIGO 
sensitivity measuring a non-resonant, NS/NS inspiral 
signal seen directly overhead with unknown jet width 
leads to a la constraint of |#o — &o/Ho\ < 2400 km. 
A four-detector network would improve this limit to 
|0 O - #0 /-Ho I < 1200 km. A BH-NS merger, seen at com- 
parable sensitivity over a longer baseline, would reduce 
the limit to \6 - 9 Q /H \ < 300 km. All of this has as- 
sumed an optimal search strategy, but with a suboptimal 
one, i.e., using pure-GR templates, one would only reach 
a bound of \$ - 9 /H \ < 600km for BH/NS events or 
< 2400 km for NS/NS ones. 

Since the maximum (9o,9q) allowed by these bounds 
would lead to a resonance in band, the absence of such 
a resonance could allow (9q,9q) to be constrained much 
better. In fact, this could constrain \9q — 9/Hq\ < 50 km 
to la assuming an optimal search strategy (see e.g., Fig. [2] 
for 2<t constraints). If wc additionally assume a four- 
detector coincident GW/GRB measurement of a NS/BH 
signal, the absence of a resonance could lead to a bound 
of 3 \9o — 9o/Hq\ < 25 km, even allowing for an unknown 
jet width and a suboptimal search strategy. 



B. Population constraints 

Combining TVo GW/GRB observations will improve 
the CS threshold roughly as y/No. Birefringence would 
manifest as a bimodal distribution of pj p, with peaks on 
either side of unity. If GRB jets are poorly collimated, 
both peaks of the distribution of p/p would be smeared 
downward due to the distribution of £ (see above) . 

Advanced detectors will observe very many short 
GRBs with GW associations. Roughly R s k y — 100 — 200 
short GRBs occur on the sky per year that produce 
enough flux to be detected by the current generation 
of satellites |59l - l61| . As a first approximation, host as- 
sociations suggest that the satellite-detected short GRB 
population is roug hly linearly distributed in redshift out 
to z ~ 0(1) [5a, [SE [fLU, [S3- Therefore, adopting an 



3 This constraint is an optimistic one, since if this bound is satu- 
rated, the resonance occurs near 2 kHz. At such high frequencies 
our calculations is limited by systematic errors due to incorrect 
waveform modeling, since one should in principle include the 
plunge and merger parts of the waveform in the modeling as 
well. 
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isotropic network horizon 4 distance D gw > 200Mpc, 
large enough to enclose many binary mergers per year, 
roughly 
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satellite-detected short GRBs could be detected by a GW 
interferometer network. Even under more conservative 
assumptions including the local event rate and expected 
beaming angles, a significant number of NS-NS or BH- 
NS binaries should be coincident with a short GRB event: 
more than > 0(10%) of all GW detections. Given the 
latest event rate predictions [65|, [66|, at least 2 should 
be detected per year. To summarize, a population of at 
least 4 short GRBs could improve constraints on 9q by a 
factor \fi = 2. 

Combining this factor of two and those of the previous 
subsection, we find that parity violation could be con- 
strained to the level of (v) < l/(p^/N d N ) for a NS/BH 
event. Modeling parity violation with CS gravity and 
first adopting a non-resonant parity-violating event, it is 
reasonable to expect that after two calendar years (the 
length of the LIGO S5 run) an advanced LIGO- Virgo 
network operating in coincidence with GRB observato- 
ries would yield a ler confidence detection of or upper 
limit on \9 - 9 o /H \ of 150 km for BH-NS (600 km for 
NS-NS) using optimally matched templates, or 300 km 
for BH-NS (1200 km for NS-NS) using a conventional 
pure-GR data analysis strategy. Again, as the maximum 
(0q,9q) allowed by the aforementioned constraints would 
lead to a resonance in band, the absence of resonances 
could allow a much better constraint. In the neighbor- 
hood of 6*o = 0, the absence of such resonance could yield 
a bound on 9q of O(10 km). 

The projected constraints on CS gravity can be thus 
more than an order of magnitude better than current 
Solar System bounds placed with the LAGEOS satel- 
lites [38[ . Although none of these constraints are compet- 
itive with the one recently placed with the double binary 
pulsar |39| > the latter only samples parity violation in the 
neighborhood of this binary, while the Solar System one 
is sensitive to parity violation only in the vicinity of the 
Earth-Moon system. The test proposed here would allow 
for constraints of the dynamical sector of parity violating 
theories along the GW geodesic, thus exploring a much 
larger region of the Universe. 



4 Because each short GRB is presumed optimally oriented, the 
range to which it can be detected is significantly greater than 
the angle-averaged range to which a typical source could be ob- 
served. In practice, the distance to which short GRBs could be 
detected will be greater still, because when a shorter stretch of 
data is searched, the detection threshold can be set significantly 
lower [Hil. 



C. Temporal-Spatial Variability of Parity Violation 

In this paper, we have proposed limits on 9q and 9q as- 
suming (i) that 9 varies on cosmological timescales and 
(ii) homogeneously and isotropically throughout the Uni- 
verse. In principle, neither assumption need apply, as the 
9 field that parametrizes parity violation in CS gravity 
should be prescribed by the more fundamental theory of 
which CS gravity is only an effective model. If the 9 field 
does not obey the symmetries of the background (homo- 
geneity and isotropy in the cosmological case), then the 
background metrics employed, on which GW propagate, 
need not be solutions to the effective theory. More gen- 
eral 9 fields would possibly require a re-analysis of the 
appropriate cosmological solutions and the propagation 
of GWs in them. 

With this caveat in mind, the more generic expressions 
of Eq. (|29|) that parameterize generic parity violating ef- 
fects via the amplification process of Sec. |V] suggest that 
the proposed constraints are primarily on a ratio of .Fr,l, 
with one factor at the source (the numerator) and an- 
other, more critical one at the detector (the denomina- 
tor). These two factors can be separated. For exam- 
ple, the resonant frequency and resonant bandwidth in 
our detectors are set primarily by the present-day evolu- 
tion of 9. If a resonance is indeed in band, every source 
with some amount of resonant polarization will be am- 
plified roughly the same way; stacking power from mul- 
tiple detections could constrain the presence of a univer- 
sal spectral feature and potentially directly constrain 9q, 
independent of its cosmological behavior, to lie outside 
the wavelength region to which LIGO is sensitive. On 
the other hand, the difference between the values of 9 
at present and the source determine how close a zero in 
the numerator is to the resonance and thus the ampli- 
fication Gr,l', cf. Eq. (p5|) . Therefore, if the detected 
SNRs do exhibit some amplifications p/p versus D that 
are not roughly uniformly increasing with distance, these 
measurements can be fit to either a functional form (if 
the points trace a curve) or a distribution (if the points 
have excess scatter) for 9. As a concrete example, so long 
as resonances are at frequencies far higher than in band 
(k9' <C 1) and the source is at low enough redshift that 
cosmological terms can be ignored (a = a s ra 1), the rel- 
ative amplification in SNR due to CS birefringence is [cf. 
Eq. (ED] 



p/p « l + A R>L 7r(/)(^-0, 



(78) 



which applies to any 9' [jf) which varies on scales faster 
than cosmological timescales. This would then allow for 
generic tests of gravitational parity violation, irrespective 
of the specific theoretical model. 
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The emergence of secular terms in the solution to cer- 
tain differential equations signals the breaking of uniform 
perturbation theory. In order to exemplify this concept 
in more detail, consider Duffing's differential equation 



ij + y + ejf = 0, 



(A.4) 



with e < 1 and the boundary conditions [y(0), y(0)] = 
[1,0]. Although the e y 3 term cannot be interpreted as 
a driving force with a certain natural frequency, we shall 
show next that its perturbative solution unavoidably in- 
troduces secular growth, just as in the case of the driven 
harmonic oscillator at resonance. 

Let us assume the following perturbative ansatz for the 
solution to Duffing's equation 



Appendix: Singular Behavior and Resonances in CS 
Gravitational Waves 



Parity violating GWs need not present a pole in the 
frequency domain, but the theory we use to exemplify 
such violation, CS gravity, does. As shown in Sec. [TT1 
such singularities arise due to poles in the modified field 
equations, i.e., Eqs. (JTUJ) and ([TBI . The effect of such 
singularities was analyzed in Sec. |V] where it was seen to 
lead to a resonant magnification or amplification of the 
SNR. That section, however, does not explain why such 
amplification is finite; this is the topic of this appendix. 
We begin with a basic introduction to resonances and 
the non-uniform validity of perturbation theory, as this 
topic may be unfamiliar to the GW community, following 
mostly [67[ • We then proceed to the study of singularities 
in CS gravity, their characterization and treatment. 



1. Resonances and the Nonuniform Validity of 
Perturbative Expansions 

Consider the differential equation for a driven har- 
monic oscillator 



y + y = cos uit , 
whose solution is 

y{t) = A cos t + B sin t + 



coswt 

1-CJ 2 



(A.l) 



(A.2) 



This solution is valid only when |c*i| 7^ 1, as in the limit 
u) — V 1, the oscillator absorbs large quantities of energy 
and the amplitude of the oscillations grow unbounded. 
The solution to Eq. (|A.1|) at resonance (|w| = 1) is 
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A cos t + B sin t 
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■ sint 
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(A.3) 



Notice that the divergent term has been replaced by a 
regular, secular contribution. Such a secular term grows 
with time and is unbounded in the limit t — ¥ 00, as then 
the oscillator absorbs energy without limit. 



y(t) = Xy v „(t). 



(A.5) 



Order by order in multiple scale analysis, Duffing's equa- 
tion becomes j/o + 2/o = to O(e ), and jji+yi = —y^ to 
0(e). The solution to the zeroth-order equation that sat- 
isfies the above initial conditions is simply yo(t) = cost. 
The solution to the first-order equation that satisfies the 
same initial conditions leads to the full solution 



y(t) = cost + 



3t 



(cos3t-cost) sint e + 0(e 2 ). 

(A.6) 

Notice that a secular term has appeared in the first-order 
solution, i.e., the last term inside the square-brackets is 
unbounded as t — > 00. This solution is valid to C(e 2 ) for 
fixed t, but as one considers times t = 0(l/e) or larger, 
Eq. (|A.6[) ceases to be valid. 

The non-uniform validity of the perturbative solution 
leads to an inaccurate approximation to the exact so- 
lution at late-times, as the exact solution is actually 
bounded for all times. To show this, one can construct 
an energy integral and show that it is bounded. In the 
case of Duffing's equation, one can multiply Eq. (|A.4[) by 
y and then rewrite it as 



1 



1 



;V 



; e v 



C, 



(A.7) 



after one time integration with constant C . Equa- 
tion (|A.7|) can be thought of as a closed bounded orbit 
in the phase space with coordinates (y,y). Due to the 
initial conditions, one easily finds that C = (2 + e)/4. 
Since each term on the left-hand side is positive if e > 0, 
then y 2 < 2C and y is bounded by y < (1 + e/2) 1 / 2 . 
We have thus proved that although the perturbative so- 
lution presents divergent behavior at late times due to 
the secular term arising at resonance, the exact solution 
is properly bounded. Another way to show this is to 
partially sum all the divergent, secular terms in the per- 
turbative solution. In the case of Duffing's equation, one 
would find that these terms sum up to cos[t(l + 3e/8)], 
which is properly bounded at t — > 00. 
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This subsection has taught us several important 
lessons. Resonant solutions to differential equations can 
be studied close to their singularities via perturbation 
theory. Such a study might reveal secular growth in the 
perturbative solutions, but this does not mean that the 
solution, or the theory from which this derives, is un- 
stable. Instead, such unbounded growth is an indication 
that the perturbative solution has a non-uniform region 
of validity. The construction of positive energy integrals 
allows us to prove that solutions are bounded from above, 
and that secular growth is artificial. 



2. Boundedness of CS Resonances 

Let us now consider the resonant structure of the CS 
differential equations for wave propagation and their per- 
turbative solution. In Scc.|TI]we found that the differen- 
tial equations for CS phase evolution have a singular driv- 
ing force, i.e., the right-hand sides of Eqs. (|TU|) and (fT"3")) 
diverge in the limit A n L fc# — > 1. Later in Sec. Mil we 
found that the exact solution to these differential equa- 
tions are also singular in this limit, i.e., Eq. ([25]) diverges 
as <S RjL — > 0. 

The perturbative, approximate solutions to the singu- 
lar CS phase evolution equations are in fact regular, al- 
though they lead to secular growth. This can be seen 
from Eqs. (|T6|) and (fT7|) . which are proportional to the 
GW travel time, or equivalcntly distance, so that as 
t — > oo, the perturbative solution is unbounded. Such 
secular growth is nothing but an indication of the non- 
uniform validity of the perturbative solution. In fact, 
from the structure of Eqs. (|16p and (|T51) . we can de- 
duce that the approximation seizes to be valid when 



the GW has travelled a distance D == O(f 1 8 ) or 
D = Oif^H^e, 1 ). 

Let us now try to show boundedness of the solution 
through the construction of an energy integral. Equa- 
tion (pM|) can be rewritten as 

h'^ + 2-^(lnS Rth )h' K + K 2 h K = 0. (A.8) 

For simplicity, we take h to be real and multiply the 
above equation by h' . Integrating once we find 



(/<;, 



(Kf 



^(ln5 R , L ) 



dr, + K Z hi = 2C . (A.9) 



Choosing initial conditions h K (0) = 1 and h' K = 0, we 
find that C = k 2 /2 > 0. It then follows that the solution 
is bounded from above by 



^J W 



d n c 

dr\ 



di) 



h K <Jl 

provided the modified scale factor is positive: 
s > c 1 2H-\ R ^k(6 + H6 



(A.10) 



UR,L 



The positivity of the effective Hubble parameter 
5 R . L /<S R ,l depends on the size of 8. This quantity is in- 
deed positive when 9 = = 8 and the universe is expand- 
ing such that H > 0. Moreover, this quantity is also 
positive away from the resonance, i.e., if 8 <C A RL /fc, 
as then the CS modification is a small deformation of 
the GR result. At resonance, (8 = A RL /fc, 8 = 0), the 
positivity of 5 RL /5 RL depends on whether 8 approaches 
A RlL /fe from below or above. If approached from below, 
= A R . L (1 - |e|)/fc as e -> 0+, then S R , L /«S R , L > 0, 
while the opposite is true if approached from above, 
(? = A R , L (l + |e|)//case->0 + . 

For strong CS parity violation, i.e., 8 > C(l/fc), 
<S R , L /<S R , L > for all 8 > 2A RL /fc, but we cannot prove 
the solution is bounded via this energy integral method if 
A RX fc# is in the interval [1, 2]. A full treatment of strong 
CS parity violation, however, requires that one includes 
higher-order ^-interactions in the action, as Eq. ([3]) is 
an effective theory valid in the weak coupling limit only. 
Such interactions will contribute to the energy integral, 
particularly in the strong coupling limit. Since we do not 
have a higher-order completion of CS gravity, we cannot 
prove boundedness in this limit. 

In spite of the apparent secular growth of the perturba- 
tive solution in the weak CS coupling limit, we have here 
shown that the exact solution remains bounded in this 
limit. The secular growth observed is then nothing but 
an artifact of the perturbative expansion and should not 
be seen as an indication of an instability in the theory. 



3. Trans-Singular CS Resonances 

Proving that the exact solution to the CS phase evo- 
lution equation remains bounded does not tell us how 
to deal with its singularity. However, as in the case of 
the driven harmonic oscillator of Eq. (|A.1[) . the solution 
we found away from resonance is not valid at resonance. 
In order to find the solution close to the singularity, we 
must solve the phase evolution equation perturbatively, 
restricting attention to a neighborhood close to the pole. 

The general form of a parity violating wave equation 



"r.L.K + K "R,L,K + 



F 

a, 
F. 



■h' = 



(A.12) 



has a pole when ,F R>L vanishes. Equation (|A.12[) is noth- 
ing but Eq. (|22|) with the replacement % — > <S R L /<S RL = 
•^r l/C^-^r.l), where wc recall that the squared amplitude 
,F R>L = <S R L . With this notation, recall that in CS gravity 
J"r.l = a" — A RL K(9'. Let us linearize the source function 
.Fptx about the (ft-dependent) conformal pole rj p (n), de- 
fined via a 2 (r]p) — K0'(r) p ) = 0, though the Taylor expan- 
sion .Frl = J-'nxlp + J'iiApiv ~ Vp)^ where the vertical 
bar stands for evaluation at the pole. Equation (|A.12() 
simplifies to a Bessel equation 



<5r,i 



A R , L fc0 



>0. (A.ll) 



K,l, k + (v- v P ) l K + K 2 /lR, L ,K = o. 



(A.13) 
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Notice that by linearizing about 77 = rjp we have lost 
the distinction between left and right-polarized waves, as 
the A R-L dependence in the source term appears only at 
second order in the Taylor expansion. 

Solutions to this wave equation linearized about the 
pole can be straightforwardly obtained as a superposition 
of Hankel functions. Imposing a no-incoming-radiation 
boundary condition at infinity, the solution reduces to 
the outward propagating, cylindrical wave solution 



ft.R,L,«(»?) = C-Ho [k(v-Vp)] 



(A.14) 



where C is a normalization constant. The quantity Hq 
is a Hankel function of the first kind, not to be con- 
fused with the Hubble parameter at present time. No- 
tice that this solution is discontinuous and singular at 
77 = rjp, a conformal time often reached after propaga- 
tion on cosmological scales. Guided by analytic contin- 
uation, by conserved current across the singularity (the 
Wronskian of Eq. (|A.13[1 ). and by the need for a nearly- 
Minkowski gravitational propagator in the neighborhood 
of r\ ps rjp(f), we adopt this particular solution on both 
sides of the resonance. 

One can reinterpret the solution in terms of propa- 
gators. In the neighborhood of some r) s , e.g., near the 
source, each spatial Fourier component is evolved for- 
ward with a simple scalar propagator: 

hn.-LAv) = A'k(?7 i 7 7s)^rx,k(?7 s ), (A. 15) 

K^w.) = ^^-^l (A.16) 

#0 i K (Vs -Tip)} 

The functional form of this equation resembles the so- 
lution to the differential equations of para-axial optics, 
where as initial data we require only the initial state 
^rx.k( 7 ?s)- In fact, except for a fraction of a wavelength 
in the neighborhood of the pole r) ~ r\ p , the propagator 
can be well-approximated at late times by 



V v-v P 

(A.17) 
which can be obtained by either expanding the Hankel 
functions in k(tj — rj p ) 3> 1 or solving Eq. (|A.13[) in this 
limit. 

The functional form of the propagator in terms of the 
source function J-r L , namely 



V FnAv) 



(A.18) 



was already derived in Sec. IIIII Taylor expanding the 
source function about resonance 

■7"r,l(£) ~ l — Xs.,Lkpd — Xs,,j,kpd(t — tp) = —\ R , h k$(t — t p ), 

(A.19) 
where in the second equality we have used the resonance 
condition 8 = X RL /k, we find that Eq. (JA.18J) reduces 



exactly to Eq. (|A.17|) . Notice however that Eq. (|A.18|) 
is formally valid away from the pole, while the true rep- 
resentation of the solution close to the pole is given by 
Eq. (TA~T6l) . 

All of these solutions diverge when n = r) p , but the 
width of the resonance is extremely narrow in frequency 
space. This width can be approximated as 5f p ~ fpdt p , 
which for a minimal- width wave packet with St p 5f p ~ 1, 

'1/2 

implies Sf p = 0(f p ). Using that f p = 1/(2x6) and 
f p = -6/(2n8 2 ), we find that 5f p ps 0[{H/6) l/2 l where 
we have assumed 9 evolves on cosmological scales 6 = 
O(H0). Evaluating these expressions today on Earth, 
we find that 6f p = O(10 _8 Hz), where for 6 we have sat- 
urated the Solar System constraint 9 < 2000 km. 

Although any single, large-scale-coherent Fourier mode 
is formally divergent, any finite-duration, band- limited 
waveform always remains finite. This is because to obtain 
any finite-duration waveform one must undo the Fourier 
decompositions performed in this paper [see e.g., Eq. ((SJ) 
or Eq. (fTTj)] . The reconstruction of the full waveform then 
requires the integration of h K (rj) over all wave-numbers: 



KAri,x l ) 



1 



(2tt) 3 / 2 
V 

(2tt) 3 / 2 



dn K^ K (r))e- iKiX , (A.20) 



where I? is a new constant related to A in Eq. (|34[) , which 
depends on the chirp mass, inclination angle and distance 
to the source, and where we have assumed the wave prop- 
agates in the x direction with conformal wavenumber k. 
This integral can be split into three integrals: two pieces 
that avoid the singularity of the propagator close to the 
pole and one piece that integrates through this singu- 
larity. Clearly, the first two pieces are finite, while the 
last piece can be approximated using Eq. (JA.16I) in the 
neighborhood of the pole: 



^p,R,] 



V 



(2tt) 3 / 2 



H (i) 

dK ° 



K(r] -rip)]e- iK ^-^ 



Hq [k (r] s - n p ) 



.7/6 ' 
(A.21) 



where we have defined k J = n p =p 2ire. An asymptotic 
expansion of the Hankel function about zero argument 
reveals that 



Hf'(x) ~ sign(x) + 



2/ 



In 



7e 



x< 1 



(A.22) 
where 7e = 0.5772157 ... is the Euler constant. Using 
this expansion, assuming the light-cone condition \ = Vs > 
and pulling out of the integral all quantities that remain 
roughly constant near the pole, the integral of Eq. (|A.21|) 
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becomes 

"■p,R,L 



h',p 

2/ 

7T 



dn {sign(Kp(?7o - r) p )) 



In 



£ = 



V 



I«p Oto - r? p ) 



-1/6 



•7E 



(2*)"* itf^ (•?.-%)] 



(A.23) 



Using ?7o = t] p +1Z(k— k p ), wc find that the most divergent 
term of the above integral is proportional to 

/ip,R, L ex £ — In ( Kp TZe) = £ — In (ff) , (A.24) 

Kp Kp yOKpJ 

where in the second relation we have used that 1Z = 
0(6k~ 2 ). We have then proved that any finite-duration, 
band-limited waveform always remains finite, as h PtRL 
remains finite when e — > 0. Furthermore, the contri- 
bution of the inverse Fourier transform near the pole is 
clearly subdominant: to use a numerical example, satu- 
rating the Solar System constraint 9q = 2000 km, such 
that k p ~ 150 Hz, and using a buffer of size e = 10 _8 Hz 
around the pole, then h p < 10 _8 £ . 

One might worry that although the pole does not af- 
fect the waveform itself, it might affect the calculation of 
the SNR, as this depends on the integral of the square 
of the Fourier transform of the waveform. The approxi- 
mate temporal-only Fourier transform of the propagator 
of Eq. (|A.16[) in the stationary phase approximation is 



K K {u\ m ) = ch 1 K K (ri,fio)e ,ur i 



I Hm) 
Hvs) 



S(LU~K)e l ^° , 



(A.25) 
where uj is conformal frequency. Then, the present-day 
temporal-only Fourier transform h(uj) can be expressed 
in terms of the spatial Fourier transform h K at the fiducial 
source time tjq: 



h(u) = 



pK K (u;) h K (r, ) e- iK * , (A.26) 



where rjQ—rj s = x- Except for an extremely narrow neigh- 
borhood near tjq ^ r)p(f), the stationary phase propaga- 
tor and thus this model for h(u>) is an excellent approx- 
imation. Substituting this expression into Eq. (|30[) for 
SNR gives the expression for resonant SNR adopted in 
the text [Eq. (|50"])]. 

As the denominator in the SNR inevitably includes a 
pole at the zero of T K (k = A RL /#'), this approximate 
integral indeed formally diverges logarithmically. How- 
ever, this logarithmic divergence arises by applying an 
approximate propagator and its Fourier transform out- 
side of their domain of applicability. In the neighbor- 
hood of /p. s , the relevant propagator is approximately 
Eq. (|A.16|) . modulo higher-derivative corrections to T' IT 
[cf. Eqs. (IAT21) . (|A~T3l) ]. 



The temporal-only Fourier transform of the resonant 
propagator is proportional to the Fourier transform of the 
standard Hankcl function, keeping in mind this Fourier 
transform (unlike many used in classical scattering prob- 
lems) is over a two-sided argument: 



K(w) 



a iu)r\ 



nPMv-rb)] 



H o [ K iVs-V P )} 



drj , 



« H^( K (r, s -r, p )) 

/OC 
dXH^(X) 
-oo 



XY 



The numerator, the Fourier transform H, 



(i) 



(A.27) 



(A.28) 



is well- 



behaved and unaffected by the proximity of a; or k to 
resonance (770 = Vp( K ))- The resonance condition en- 
ters only in the denominator. As Hq (X) has no zeros 
for real X, however, the time-domain Fourier transform 
of the resonant propagator K and the transform of the 
amplified waveform h(ui) have no poles on the real axis. 
This shows that even if the calculation where done ex- 
actly, the true SNR would remain finite. Moreover, given 
the extremely slow rate at which the resonant frequency 
evolves, we can approximate Hq (u)/ k) as proportional 
to 6(ui — k). This implies that to a good approximation 
the Fourier transform of the time dependent propagator 
is K K (r]o,r] s )d{u! - k). 

With this reasonable approximations at hand, the SNR 
becomes 



Pr,l 



U<-\ 



|/Ir,l 



v- 



Sh(n) 
dn 



-dK , 



\K K ( Vo ,Vs)\ 2 K-y 3 - 

bh(K) 



(A.29) 



As before, let us split this integral into three pieces: two 
of them that avoid the singularity and one that integrates 
through it. In the neighborhood of the singular point, we 
once more approximate the propagator using Eq. (|A.16I) , 
so that 



R,LjP 



4£H 



r(l) 



H o [k{vo-V p )} 



H o [k(V3 -rip)] 



.-7/3 _ 



Sh(n) ' 
(A.30) 



Notice that the argument of the Hankcl function func- 
tion in the numerator is the only quantity that formally 
diverges at the pole, so let us now rewrite the integral as 



2 

rR,L,p 



£ 



£ 



4V 2 



r(i) 



#6 [KpiVo ~ Vp)] 



H o [ k p(Vs -V P )] 



dK , 

.-7/3 
Kp 

Sh(Kp) 



(A.31) 



where we have pulled out several factors that are non- 
singular at the pole and roughly constant near the res- 
onance. These expressions lead to singular results if 
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one evaluates it with the late-time expansion of the 
Hankel function of Eq. (|A.17|) , Instead, using the 
small- argument Hankel function expansion provided in 
Eq. (|A.22j) . the SNR integral becomes 



(A.32) 



Pr, l ,p - £ / {sign[K p ^(K - k p )] 



[ ' KpU ^ Kp)l | 



d,K . 



whose most divergent contribution is asymptotic to 



pL, p ~ £ e [ m (KpTle)} = £ e 



In 



8nl 



(A.33) 



Notice again that the contribution to the SNR near the 



pole is finite and in fact it formally vanishes in the limit 
8n p — > 0. As is clear from the above analysis, regulariza- 
tion is here justified: for a suitably small rcgularization 
interval e, the contribution from a neighborhood e around 
the pole is small. To provide a numerical example, satu- 
rating the Solar System constraint again (k p = 150 Hz) 
and choosing a regularization interval of size e = 10 _8 Hz 
we find that p R-LiP < 30/5 RLiP , where /S RL ,p is the con- 
tribution of the SNR in an interval of size 5n p in the 
absence of a resonance. The SNR is mostly dominated 
by contributions between 100-300 Hz, so splitting this 
region into segments of size 6f p = 10~ 8 Hz, one obtains 
approximately 10 10 segments. Assuming each segments 
carries approximately the same power, one can estimate 
that /5 R , L , P - 10/10 10 = O(10" 9 ). Therefore, the contri- 
bution from the resonance that we lose by regularizing is 
approximately p RL .p ^ 10 -8 and can be neglected. 
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